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The phase structure of Wilson fermions with twisted mass (/i) has been investigated in [JTJ, 
As it is explained there, the observed first order phase transition limits the minimal pion mass 
which can be reached in simulations at a given lattice spacing: m™ n ~ G(a). The phase structure 
is schematically depicted in the left panel of Fig. [Ij The phase transition can be observed in 
simulations with twisted mass fermions, for instance, as a "jump" or even metastabilities in the 
average plaquette value as a function of the hopping parameter (fc). (The right panel of Fig. 
shows thermal cycles to demonstrate this phenomenon.) 

One possibility to weaken the phase transition and therefore allow for lighter pion masses at 
a given lattice spacing is to use an improved gauge action like the DBW2, Iwasaki, or tree-level 
Symanzik (tlSym) improved gauge action instead of the simple Wilson gauge action. This has been 
successfully demonstrated in [Q, [|, ||] . 

Here we report on our attempts to use a smeared gauge field in the fermion lattice Dirac 
operator to further reduce the strength of the phase transition. This is relevant in simulations with 
Nf = 2 + 1 + 1 (u,d,s,c) quark flavours ^ where the first order phase transition becomes stronger 
compared to Nf = 2 simulations. The main impact of the above mentioned improved gauge actions 
on the gauge fields occuring in simulations is to suppress short range fluctuations ("dislocations") 
and the associtated "exceptionally small" eigenvalues of the fermion matrix. The same effect is 
expected from smearing the gauge field links in the fermion action. The cumulated effect of the 
improved gauge action and smeared links should allow for a smaller pion mass at a given lattice 
spacing and volume. Our choice is the Stout smearing procedure as introduced in since it 
can easily be implemented in the Hybrid Monte Carlo (HMC) based updating algorithms we are 
currently using. 

One should keep in mind that a possible caveat of this procedure is "oversmearing", i.e., re- 
moving too many small eigenvalues by applying too many smearing steps and/or using a too high 
value for the smearing parameter — because not every small eigenvalue is "unphysical". In addi- 
tion, after many smearing steps the fermion action can become too delocalised which can lead to 
an unwanted slowing down of the approach to the continuum limit. In order to avoid this caveat we 
choose to work with only one step of very mild Stout smearing. Moreover we keep these smearing 
parameters fixed as we change the lattice spacing. 

In Section [I] we will shortly review the smearing procedure and the twisted mass formulation, 
as well as some details concerning the used updating algorithms. Section || is devoted to the 
presentation of the results of our numerical simulations using Nf = 2 and Nf = 2 + 1 + 1 flavours 
of twisted mass quarks. 

1. Stout smearing and twisted mass fermions 

1.1 Analytic smearing for SU(3) link variables 

To have a smearing procedure which is analytic in the unsmeared link variables is an essential 
feature, if one is to use the smeared variables in an updating scheme like HMC based algorithms 
requiring the calculation of the derivative (or force) with respect to the unsmeared link variables. 
The Stout smearing procedure as introduced by Morningstar and Peardon in [|7|] was designed to 
meet this requirement. We will briefly describe it in the following but will limit ourselves to the 
case of SU(3). For more details we refer to the cited work. 
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The (n + l) th level of Stout smeared gauge links is obtained iteratively from the n th level by 



u r> u . 



(i.i) 



In the following we will refer to the unsmeared ("thin") gauge field as = uj?\ while the maxi- 
mally smeared ("thick") gauge field for Af-level Stout smearing will be denoted as = Uu ■ The 
SU(3) matrices are defined via the staples C M : 



Qt\x)= { -W{x) C f{x)-h.c. 



+Ui n) \x - v)U^ (x - v)ul n) (x - V + ft) ) , 



r(n), 



:Tr 



U^\x)cf{x)-h.c. 



(1.2) 



(1.3) 



where in general p^ v is the smearing matrix. In our numerical simulations we used exclusively 
isotropic 4-dimensional smearing, i.e., p^ v = P- 



The thick gauge field will only be used in the fermion operator, cf. Eq. (|1.7|). Therefore the 
usage of smeared links is nothing else but a different discretization of the covariant derivative 
operator on the lattice. For the gauge part of the action the thin gauge field still is the relevant one. 



1.2 Twisted mass fermion action 

The notations in this subsection follow Ref. [^]. (For details we refer to this work.) We 
performed simulations with one light doublet (u, d) of twisted mass Wilson fermions (only using 
the fermion matrix Q) ). Later on we shall add a second doublet for the heavier quarks (c,s), 
where the masses are non-degenerate due to the addition of an extra mass term (cf. ||8|, |9]]). The 
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fermion action then reads 



nfermion 



qI x) 



■V-.V 



yCixQucyXhy + %h,xQ*xy% h ,y) ' 



1 



±4 



(1.4) 

(1.5) 
(1.6) 

(1.7) 



where pL Kx = l/(2Kx) is the untwisted mass and a/i/, a/i/j, and apLg are the twisted mass terms in 
the light and heavy doublet and the split mass term, respectively. 

1.3 Algorithms 

We used two different algorithms with independent implementations of the Stout smearing 
routines to be able to cross-check our results. The first algorithm is the HMC algorithm with 



multiple time scale integration and mass preconditioning as described in [10]. In that case, the 



smearing routines were taken from the CHROMA-code package [11] and a chronological inverter 
was included, too. 

Since this algorithm only allows to simulate an even number of fermion flavours (also exclud- 
ing the case of a split doublet as described in Sec. 12), as a preparation for the Nf ■ = 2 + 1 + 1 



simulations, we also added Stout smearing routines to our existing Polynomial HMC (PHMC) 



[12, 13] update code, where we perform one stochastic correction step at the end of a trajectory. 



For details on the implementation of the PHMC, cf. [ 14, 15]. We used trajectory lengths of 2 x 0.35 
to 3 x 0.35 and determinant breakup of hb = 2. 



2. Numerical simulations 



2.1 Nf = 2 

In all of the simulations presented here, we used the tree-level Symanzik improved gauge 
action on either 12 3 x 24, 16 3 x 32, or 24 3 x 48 lattices. We compare results obtained using one 
level of Stout smearing (N = 1) with p = 0.1 or p = 0.125 to simulations without smearing of the 
link variables in the Wilson twisted mass fermion action. Our choice of mild smearing (smearing 
only once with a small parameter) should guarantee that the fermion action remains well localized 
on physical scales even on relatively coarse lattice spacings. 

Figures ||| and |3] show the average value of the (thin) plaquette without and with smearing, 
respectively. In the case without smearing a jump in the average plaquette value is clearly visible. 
Here we also observed metastabilities, which show up as differences between runs starting from a 
random (hot) or ordered (cold) configuration (red circles and blue triangles, respectively, in Fig. |2|). 
In the case of Stout smearing it is unclear whether there is still a phase transition at all, since the 
left panel of Fig. || shows a rather smooth dependence of the average plaquette value on the inverse 
hopping parameter. To examine if metastabilities may still arise with Stout smearing, we started 
runs from either a random (hot) or ordered (cold) configuration at the same parameters where the 
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Figure 2: Average plaquette value without Stout smearing with hot (red circles) and cold starts (blue 
triangles) at two different values for the gauge coupling j3 and twisted mass ajj.. 
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Figure 3: Simulations with Stout smearing. Left panel: Average plaquette value on two different lattice 
sizes. Right panel: Monte Carlo history of the plaquette value for hot (red) and cold (blue) starts. 



hopping parameter was chosen to lie in the region of fastest increase of the average plaquette: 
K = 0. 1513 or 1/ (2k) 3.305; the Monte Carlo histories of the two runs are displayed in the right 
panel of Fig. [| One can see that after roughly 500 trajectories both runs thermalized at the same 
average value for the plaquette giving no evidence to the existence of metastabilities. Figure |] 
shows the (untwisted) PCAC quark mass (left panel) and the squared pion mass (right panel) as 
a function of fj. K . The former also shows no clear evidence for the presence of a phase transition, 
since both branches (positive and negative PCAC quark mass) extrapolate to roughly the same 
critical value of pL K . From the latter one can read off that on the volume L/ro ~ 4 a minimal pion 
mass of m K ro ~ 0.7 is easily achieved for a ks vq/4. (We use here for setting the scale the Sommer 
parameter ro.) 

On the 24 3 x 48 lattice in a previous simulation at /3 = 3.8 [16] without smearing we observed 
a problematic behaviour in time histories implying very long autocorrelations. We are presently 
repeating this run with one level of Stout smearing (p = 0.125) to see the effect of smearing for 
such a situation. Although our first results indicate that smearing helps, it is too early to give a 
definite conclusion at this point. It will also be interesting to compare physical observables, e.g., 
fps and mps, between the Stout smeared and unsmeared simulations. 
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Figure 4: Left panel: PCAC quark mass, right panel: squared pion mass using Stout smearing 



© 



0.52 



0.51 



0.5 - 



© 



N f =2+1+1, 12 J x24 lattice 
tISym gauge action 
B = 3.25 
an, = 0.01 
au„ = 0.315 
an 6 = 0.285 



®2© 



© 



© 



2.86 2.87 
1/(2k) = 1/(2k,) = 1/(2^) 



Figure 5: Average plaquette with Nf = 2+ 1 • 
j3 = 3.35 (right panel) lattices, see [g]. 
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2.2^ = 2 + 1 + 1 

Recently, the possibility of adding the strange quark in dynamical twisted mass simulations 
has been tried following the lines of [^] by introducing a mass splitting term in the heavier doublet, 
see Eq. (|L6|). In that way not only a strange quark will be added but also the much heavier charm 
quark is taken into account. For first numerical results see ^J, where an important conclusion is 
that the extra dynamical quarks strengthen the first order phase transition. As an example, in Fig. || 
we show the jump in the average plaquette for two different lattice spacings at a fixed physical 
volume. On the coarser lattice spacing (left panel) again metastabilities show up. At the finer lattice 
spacing (right panel) there are no more metastabilities but there is still a considerable "jump" in 
the average plaquette. The findings from our Stout smeared run for Nf = 2 suggest that smearing 
could substantially help in the case of Nf = 2 + 1 + 1. 



Conclusions & Outlook 



The conclusion of testing Stout smearing with twisted mass Wilson quarks is that the first order 
phase transition at non-zero lattice spacing becomes weaker as a result of smearing. Therefore 
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moderate Stout smearing can be an option — in particular for future numerical simulations in the 
twisted mass formalism with dynamical u-, d-, s- and c-quarks. 
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